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Periodic homogenization with an interface 
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We consider a diffusion process with coefficients that are periodic outside of an 
'interface region' of finite thickness. The question investigated in the articles 1 ' 2 
is the limiting long time / large scale behaviour of such a process under diffusive 
rescaling. It is clear that outside of the interface, the limiting process must 
behave like Brownian motion, with diffusion matrices given by the standard 
theory of homogenization. The interesting behaviour therefore occurs on the 
interface. Our main result is that the limiting process is a semimartingale whose 
bounded variation part is proportional to the local time spent on the interface. 
We also exhibit an explicit way of identifying its parameters in terms of the 
coefficients of the original diffusion. 

Our method of proof relies on the framework provided by Freidlin and 
Wentzell 3 for diffusion processes on a graph in order to identify the generator 
of the limiting process. 



1. Introduction 

In this note, we report on recently obtained results 1 ' 2 on the long-time 
large-scale behaviour of diffusions of the form 

dX = b(X)ds + dB(s) , X(0) = xeM. d , (1) 

where B is a ci-dimensional standard Wiener process. The drift b is assumed 
to be smooth and such that b(x + e^) = b(x) for the unit vectors with 
i — 2, . . . , d (but not for i = 1). Furthermore, we assume that there exist 
smooth vector fields b± with unit period in every direction and 77 > such 
that 

b(x) = b + (x) , x\ > 77 , b(x) = b-(x) , x\ < —r\ . (2) 

Setting X e (t) = eX(t/e 2 ), our aim is to characterise the limiting process 
X = lim e — >o X £ , if it exists. In the sequel, we denote by C the generator of 
X and by £± the generators of the diffusion processes X± given by Eq. (1) 
with b replaced by b±. The processes X± will be viewed as processes on 
the torus T d , and we denote by /i± the corresponding invariant probability 
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measures. In order to obtain a diffusive behaviour for X at large scales, we 
impose the centering condition J Jd b±(x) n±(x) = 0. 

Before stating the main result, first define the various quantities involved 
and their relevance. 

We define the 'interface' of width 77 by T v = {x e R d : x\ € [—77,77]}. 
In view of standard results from periodic homogenization, any limiting 
process for X s should behave like Brownian motion on cither side of the 
interface X = {x\ — 0}, with effective diffusion tensors given by 

D% = I (5 ik + d k gf){5 kj + d k gf) d^ ± . (3) 

(Summation of k is implied.) Here, the corrector functions g± : T d — > M. d are 
the unique solutions to C±g± = —b±, centered with respect to /x±. Since 
b± are centered with respect to /j,±, such functions do indeed exist. 

This justifies the introduction of a differential operator C on M. d defined 
in two parts by £+ on I + = {x\ > 0} and £_ on 7_ = {xi < 0} with 

D ± 

C± - -fdtdj , (4) 

then one would expect any limiting process to solve a martingale problem 
associated to C. However, the above definition of C is not complete, since 
we did not specify any boundary condition at the interface Jo- 
in the one dimensional case 1 the analysis is considerably simplified since 

• The interface is zero dimensional in the limit and hence cannot 
exhibit any more complicated behavior than preferential exit be- 
havior. 

• The non-rescaled process is time-reversible and therefore admits an 
invariant measure for which one has an explicit expression. 

Both of these clues allow us to make a reasonable guess that in one dimen- 
sion the limiting process will be some (possibly different on each side of 
zero) rescaling of skew Brownian motion. Since the diffusion coefficients on 
cither side of the interface are already determined by the theory of periodic 
homogenisation, the only parameter that remains to be determined is the 
relative probability of excursions to either side of the interface. This can be 
read off the invariant measure by using the fact that the rescaled invariant 
measure should converge to that of the limiting process. 

One of the main ingredients in the analysis of the behavior of the limiting 
process at the interface is the invariant measure /1 for the (original, not 
rescaled) process X. If we identify points that differ by integer multiples 
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of ej for j = 2, . . . , d, we can interpret X as a process with state space 
R x T d_1 . It then follows from the results in 5 that this process admits a 
er-finitc invariant measure [i on R x Y d ~ x . 

Note that the invariant measure /x is not finite and can therefore not be 
normalised in a canonical way. However, if we define the 'unit cells' by 

C+ = [j,j + l]xT d ~ 1 (5) 
CJ = [-j-l,-j]xT d - 1 (6) 

then it is possible to make sense of the quantity q± = lim^oo fj,(C^). 
Let now p± be given by 

+ q-D ii 

Unlike in the one-dimensional case, these quantities are not sufficient to 
characterise the limiting process since it is possible that it picks up a non- 
trivial drift along the interface. It turns out that this drift can be described 
by drift coefficients aj for j = 2, . . . , d given by 

aj =2(2±r + ^) [ bj(x)rtdx), (8) 

where /x is normalised in such a way that q + + 17 _ = 1. 

Given all of these ingredients, we construct an operator C as follows. 
The domain T>(£) of C consists of functions / : M. d — > R such that 

• / is continuous and its restrictions to I + , /_ , and X are smooth. 

• The partial derivatives dif are continuous for i > 2. 

• The partial derivative dif(x) has right and left limits dif\± as 
x — > Jo and these limits satisfy the gluing condition 

d 

P+di/|+-P-di/|-+5>^/ = 0. (9) 

For any / G f(xC), we then set Cf{x) = C±f{x) for a; e I±. With these 
definitions at hand, we can state the main result of the article: 

Theorem 1.1. The family of processes X 6 converges in law to the unique 
solution X to the martingale problem given by the operator L. Furthermore, 
there exist matrices M± and a vector K e R d such that this solution solves 
the SDE 

dX(t) = l {jtl < 0} M-dW{t) + l { x 1>0} M + dW(t) + K dL(t) . (10) 
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where L denotes the symmetric local time of X± at the origin and W is a 
standard d- dimensional Wiener process. The matrices M± and the vector 
K satisfy 

M±Ml = D ± , K 1 =p+-p_ , Kj = a, , (11) 
forj = {2,...,d}. 

In Figure 1, we show an example of a numerical simulation of the process 
studied in this article. The figure on the left shows the small-scale structure 
(the periodic structure of the drift is drawn as a grid). One can clearly 
see the periodic structure of the sample path, especially to the left of the 
interface. One can also see that the effective diffusivity is not necessarily 
proportional to the identity. In this case, to the left of the interface, the 
process diffuses much more easily horizontally than vertically. 

The picture to the right shows a simulation of the process at a much 
larger scale. We used a slightly different vector field for the drift in order to 
obtain a simulation that shows clearly the strong drift experienced by the 
process when it hits the interface. The remainder of this note is devoted to 
a short discussion of the proof of Theorem 1.1. 




2. Idea of proof 

As is common in the theory of homogenization, the pattern of the proof 
is as follows: one first verifies tightness, then shows that any limit point 
satisfies the martingale problem associated to C, and then finally identifies 
solutions to this martingale problem as the unique solution to Eq. (10). 
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2.1. Tightness of the rescaled processes 

We want to show that the modulus of continuity of X s is well-behaved 
uniformly in e. The only barrier to this holding can easily be shown to be 
the drift picked up by the process in the interface. In order to bound this, 
we thus need to show that the process does not spend too much time there. 

We decompose the trajectory for the process X e into excursions away 
from the interface, separated by pieces of trajectory inside the interface. We 
first show that if the process starts inside the interface, then the expected 
time spent in the interface before making a new excursion is of order e 2 . 
Then, we show that each excursion has a probability at least s/VS of being 
of length S or more. This shows that in the time interval 5 of interest, the 
process will perform at most of the order of VS/s excursions, so that the 
total time spent in the interface is indeed of the order eVS. Since the drift 
of the rescaled process is of order 1/e, we conclude that the modulus of 
continuity is of order VS everywhere. 

2.2. Identification of the limiting martingale problem 

In order to identify the martingale problem solved by the limiting process, 
it is possible to adapt a result obtained by Freidlin and Wentzell in the 
context of diffusions on graphs. 3 The main ingredients are the following. 
For S — e a with a € denote by t s the first hitting time of dXg 

by X £ . We then show that for p± and aj as in Eq. (7) and Eq. (8), the 
convergences 



take place uniformly over x € T ert . 

In order to show the first identity in (12), let Tk be the first hitting time 



One can then show that lim^oo \p\—p k + \ — using the fact that the process 
returns to any small neighborhood in 1^ before Tk with probability tending 
to 1 as k — > oo, allowing the process to forget about its initial conditions 
through a coupling argument. The values p± can then be computed in a 
way similar to the one-dimensional case. 

The main ingredient in this calculation is the fact that the invariant 
measure /x for the process X (which we can view as a recurrent process on 



(12) 




(13) 
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M x T d - 1 ) gets closer and closer to multiples of ji± away from the interface. 
This can be formalised as: 

Proposition 2.1. Let A denote a bounded measurable set and denote by 
[i the (unique up to scaling) invariant <r -finite measure of the process X. 
Denote furthermore by /x± the invariant measure of the relevant periodic 
process, normalised in such a way that jit±([fc, k + 1] x T d_1 ) = 1 for every 
k € Z. Then there exist normalisation constants q± such that, 



Urn (| M (A + fc)- g+M+ (^)| + |M^-fc)-9-M-(^)l) -0. (14) 



(Here k is an integer.) Furthermore, this convergence is exponential, and 
uniform over the set A if we restrict its diameter. 

In order to obtain an expression for the limiting values p± , one can now 
argue as follows. Considering the first component of the limiting process, it 
is reasonable to expect that it converges to a rescaling Y of skew Brownian 
motion. This is characterised by three quantities: its diffusivity coefficients 
on either side of the interface (we already know that they are given by Df^) 
and a parameter p + such that, setting p_ = 1 — p + , 



The invariant measure for Y is known to be proportional to Lebesgue mea- 
sure on either side of the interface, with proportionality constants q± = ji-. 



We can then simply solve this for p± . 

The second part of (12) is shown in two steps. With as before, we 
have the identity 



for any fixed starting point x in the interface. If k is large, then the process 
X has had plenty of time to "equilibrate" , so that it is natural to expect 
that ctj is proportional to J bj(x) fx{dx). The only question is: what should 
be the correct proportionality constant? 

In order to answer this question, let us assume for the sake of the ar- 
gument that bj = iV _1 lr_jv,iv] f° r some fixed but large value of N. (Note 
that the fact that the function bj appearing in Eq. (16) is given by the 
drift of the original diffusion is irrelevant to the argument, we could ask 
about the value of this limit for any function b that is localised around 
the interface.) We then have / bj(x) fi(dx) w 1, thanks to the normalisa- 
tion q + + g_ = 1. On the other hand, we know that the first component 




(15) 



ii 




(16) 
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of the rescaled process converges to skew Brownian motion described by 
the parameters p± and D u . Time-changing the process by a factor D n on 
either side of the origin, we can reduce ourselves to the case of standard 
skew-Brownian motion with parameters p± . Since this consists of standard 
Brownian motion excursions biased to go to either side of the origin with 
respective probabilities p± , this yields in this particular example 



where B is a standard Brownian motion. A simple calculation then shows 
that the term under the expectation is asymptotic to 2k, so that we do 
indeed recover the proportionality constant from Eq. (8). 

2.3. Uniqueness of the martingale problem 

Finally in order to show uniqueness of the martingale problem, we use 
Theorem 4.1 from 6 in conjunction with the Hille-Yosida theorem to ensure 
that the domain of the generator to our martingale problem is large enough. 
It is then possible to explicitly construct solutions to the system of SDEs 
given in (10) and to show that they solve the same martingale problem, 
thus concluding the proof. 
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